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Abstract 

Glueball masses with J < 7 are computed both for C = +1 and 
C = —1 using the string Hamiltonian derived in the framework of 
the Vacuum Correlator Method. No fitting parameters are used, and 
masses are expressed in terms of string tension a and effective value 
of Og. 

We extend the calculations done for J < 3 using the same Hamilto¬ 
nian, which provided glueball masses in good agreement with existing 
lattice data, to higher mass states. It is shown that 3 , 5 and 7 

states lie on the odderon trajectories with the intercept around or be¬ 
low 0.14. Another odderon trajectory with 3g glueballs of V-shape, 
corresponds to 11% higher masses and low intercept. These findings 
are in agreement with recent experimental data, setting limits on the 
odderon contribution to the exclusive yp reactions. 


1 Introduction 

There was a renewal of interest in glneballs recently, mostly connected to 
pomeron and odderon trajectories P-IO]- High-spin glneballs, e.g. 
were calcnlated on the lattice da. and the problem of odderon attracted 
much attention EllllEl- Search for odderon exchange at HERA Pd in the 
reactions yp —vr^p, yp —Tr'^vr^p,... has not provided any indications for an 
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existence of these processes at high energies^ and has lead to rather strong 
limits on the corresponding cross sections. This allows to put limits on the 
intercept of the odderon in the models which predict couplings of the odd- 
eron with hadrons |3]. Study of this problem allows to distinguish between 
the models and seriously question some of them. Therefore it seems to be 
necessary to make a more detailed analysis of our theoretical calculations jS] 
and to compare additional high-spin glueball masses. 

A theoretical study of glueballs in QCD was started in i-ra and is 
closely related to the problem of the pomeron, i.e. leading Regge pole, which 
determines the asymptotic behavior of scattering amplitudes at very high 
energies. It is usually assumed that the pomeron in QCD is mostly gluonic 
object and glueball resonances with vacuum quantum numbers and spins 
belong to this trajectory. Another interesting hypothetical Regge singularity 
is the ’’odderon”, which has negative signature and C parity and can be 
built out of at least 3 gluons. Most studies of the pomeron and odderon 
singularities in QCD are based on applications of the perturbation theory 

[HI. 

Our method is based on the QCD path integral formalism, where all dy¬ 
namics is encoded in held correlators - the so-called Field Correlator Method 
(FCM) [T3] (for a review see starting from that one can derive in the limit 
of small gluon correlation length A the relativistic Hamiltonian UTI, which 
effectively describes the fundamental (adjoint) string with quarks (gluons) 
at its ends. In this simple limit, A —0, Hamiltonian is local but nonlinear 
in p^,L^, which rehects complicated dynamics of relativistic rotating string 
(similar results are obtained in jlHji. 

One should stress, that nonperturbative (NP) approach started in [1^] and 
developed in jS], is based on the extrapolation of trajectory J{t) connecting 
2++, 4++, states ( for pomeron) and 3 ,5 ,7 , states (for odderon) 

to the physical region of scattering, t = <0. In doing so one assumes 

that no extra singularities appear is the J plane on the way to t = 0. In case 
of pomeron indeed two other trajectories, / and exist which intersect 
with glueball trajectory and therefore one should takes into account mixing 
between them. As a result in jHj a combined trajectory was calculated with 
realistic pomeron intercept. In this approach perturbative contributions are 
of subsidiary character and can shift the intercept by approximately 0.2. 

^Note that the usual uj,p - Regge poles with a„(0) « ap(0) « 0.5 give contributions to 
these reactions, but the corresponding cross sections decrease with energy as ^|t=o -g- 
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It is important to stress that the region t 0 as well as t ~ > 0, 

M ~ several GeV belongs to the primarily nonperturbative regime, where 
all intergluonic and interquark distances are large, hence the extrapolation 
of trajectories to t = 0 is inside the NP domain. 

In the case of odderon in jH] the intersection of the leading gluonic tra¬ 
jectory with G = —1 with qq trajectories does not take place, since p, a; 
trajectories have much larger intercepts, and the resulting odderon intercept 
in jH] obtained using the 3 glueball mass and assumed slope of {2TTaadj)~^ 
was predicted around -1.5. This value excludes possible odderon discovery in 
reactions yp —>• (vr^, /^, jHlIZI within the upper limits set by the authors. 

On the other hand, the perturbative (BFKL) approach jH] starts from 
another premises. It considers pomeron (odderon) as 2g (3g) system of 
Reggeized gluons having only perturbative gluon exchanges, which is justi- 
hed in the perturbative domain of small (< 1 GeV“^) interparticle distances, 
the situation which might be realized for large negative f, large s, and very 
small sizes of color dipoles exchanging pomeron. Therefore one encounters 
the problem of analytic continuation of pomeron singularity to the nonper¬ 
turbative domain of f = 0 (and realistic dipole sizes). 

Assuming this can be done, one obtains both pomeron and odderon in 
vicinity of J = 1 [Hj, see {20], |2I| and |3| for more discussion of 
odderon. 

While for pomeron this is reasonable (however the coincidence of pomeron 
with J = 1 is trivial in the limit of small a^), for odderon it presents a pre¬ 
diction of strong negative C - parity-odd contribution to different reactions, 
in particular of the type mentioned above jam. 

In a different approach(see jl] and refs, therein) a model nonperturba¬ 
tive picture for high-energy scattering was exploited yielding the odderon 
intercept aoddi^) = 1, which as well as the BFKL prediction is at odds with 
existing data. 

Therefore we feel it is necessary to clarify the situation with glueball tra¬ 
jectories and to this end to extend our previous calculations |H| to higher 
spin states, namely we calculate the masses of glueballs with C = +1, J = 
0++, 2++, 4++, 6++ and with C = —1, J = 1 ,2 ,3 ,5 ,7 . From 

those we calculate the odderon and pomeron trajectories and hnd the appro¬ 
priate intercepts. 

Another aim of our calculations is the comparison of newly found glueball 
masses with existing lattice data [Iliaiia 123 I21j- In our previous work jH) 
we have found a good agreement of all glueball masses with lattice data. 
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Recently new states, 4++ and 6'''+ have been computed [H |2] and we can 
compare those with our analytic results. We also compare our results with 
recent analytic calculations 1112312n]. The paper is organized as follows. In 
section 2 the Hamiltonian is given together with spin terms following [S] and 
masses for gg glueballs are calculated. 

In section 3 two possible conhgurations of 3g glueballs are dehned, a A- 
type and an F-type, and the resulting masses are calculated. In section 
4 glueball trajectories are obtained and intercepts are found and discussed. 
The concluding section 5 is devoted to the discussion of the results from 
the point of view of lattice and experiment correspondence, and to possible 
improvements. 


2 String Hamiltonian and spin corrections 

The Hamiltonian for the gg system was derived in jH] in the same way as it 
was done for the qq system I2SI. and can be written as 


H = Ho + AH,, AH, = Hsl + Hss + Ht (1) 

where Hq is the generalization of the spinless qq Hamiltonian, obtained in 
na. and given in nacHi 


Ha = 


Pr 


+ /i(t) + 


L{L + 1) 


r^[p + ‘^Io{p-ly^df3] 


+ 



^Idjdfd 2 , 1 

- - -r H— 

2 



t)dj3. 


( 2 ) 


Here fi{t) and are positive auxiliary functions which are to be found 

from the extremum condition na. Their extremal values are equal to the 
effective gluon energy (/i) and energy density of the adjoint string {v). 

Here aadj = fo'/nnd, and we shall always set ct/ = 0.18 GeV^ as found 
from meson Regge trajectories II71I23. 

To hnd the spin-averaged masses from ([Q), AHg = 0, one can use the 
WKB procedure developed in na for the qq case and having accuracy better 
than 5% for n = 0, which yields the values given in Table 1 (exact values 
available for L = 0 are given in parentheses). 
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Table 1 

Spin-averaged masses (in GeV) of gg states with L = 0, ...5 and n = 0,1, 2, af 
0.18 GeV^ 


L 

n 

0 

1 

2 

3 

4 

5 

0 

2.09 (2.01) 

2.65 

3.13 

3.53 

3.88 

4.21 

1 

3.20 (2.99) 

3.65 

4.03 

4.37 

4.67 

4.95 

2 

4.01 (3.75) 

4.40 

4.74 

5.04 

5.31 

5.56 


In (P) and in Table 1 the effect of perturbative gluon exchanges between 
gluons was neglected, and we consider it as a leading reasonable approxima¬ 
tion in obtaining glueball masses. In doing so we follow the argument given 
in [S], where it was shown that the Goulomb-like adjoint charge interaction 
is not formed between valence gluons, and moreover, following BFKL ap¬ 
proach, one can consider gluon exchange as an effectively not large which 
can be deduced from the relatively small shift A = ap(0) — 1 of the in¬ 
tercept, when terms 0{a^) are taken into account [TTl IH]. Therefore in the 
hrst approximation we neglect perturbative gluon exchanges, keeping them 
only in the spin-dependent terms AHs- This strategy is supported by the 
comparison with lattice data of the spin-averaged masses (see Table 4 of ref. 
i) and, separately, by the comparison of our spin splittings of masses with 
lattice data, to be discussed below. 

Spin-splitting terms AHg are considered in detail in jS] for L = 0,1, 2 and 
we list in Table 2 the resulting masses, in comparison with existing lattice 
data. 

The difference of Table 2 from the corresponding Table 6 of ref. [S] is 
that we £x cx/ = 0.18 GeV^ and recalculate lattice masses for this value 
of af. Moreover, we take for L = 0 in brackets as{eff) = 0.2 (for spin¬ 
splitting terms, since the corresponding spin interaction occurs at relatively 
small distances, while for the rest masses as{eff) = 0.3 was taken, as in [H]). 
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Table 2 


Comparison of calculated glueball masses (in GeV) with lattice data (cx/ = 
0.18 GeV^, tts = 0.3 = 0.2 in parentheses)) 


jPTT- 

^'^theory 

this work 

Mlat 

122] 

|231 

[24] 

0++ 

(1.61) 1.41 

1.53T0.10 

1.53T0.04 

1.52T0.13 

2++ 

(2.21) 2.30 

2.13T0.12 

2.20T0.07 

2.12T0.15 

0++* 

(2.72) 2.41 

2.38±0.25 

2.79±0.09 


2++* 

(3.13) 3.32 

2.93±0.14 

2.85±0.28 


0-+ 

2.28 

2.30±0.15 

2.11±0.24 

2.27T0.15 

0-+* 

3.35 

3.24±0.2 



2-+ 

2.70 

2.76T0.16 

3.0±0.28 

2.70±0.19 

2-+* 

3.73 

3.46T0.21 




One can see in Table 2 a good agreement of our calculated glueball masses 
with measured lattice values, especially for ground states. The radially ex¬ 
cited states (marked with an asterix) are some 5-7% higher than lattice data, 
which probably is the result of WKB approximation for the Hamiltonian 
Q; indeed as seen from Table 3 of ra. and our Table 1, second column, the 
exact mass eigenvalues for L = 0 and > 0 are 6-7% lower than those from 
WKB approximation. 

This agreement of theory with lattice data becomes even more striking, 
when one realizes that our calculation has no htting parameters at all, since 
string tension is a given scale parameter (one could compare dimensionless 
M/^/a, as it is done in Table 4 of [Hj) and ag, hxed at the characteristic 
value, Us = 0.3 describes only spin splitting of masses. 

Let us now discuss high spin states, not present in Table 2. The L = 
2,S = 2 states include 14++ states which are spread over the 
mass distance of 63 MeV due to the spin-orbit and tensor level splitting 
(see table 7 of |H1). Here perturbative and nonperturbative (Thomas term) 
spin-orbit interaction almost cancel each other. 

One can expect, that for higher states the splitting of the states will be 
even less, as it is observed experimentally for mesons and in what follows we 
shall neglect spin splitting for states with L > 2. 

The resulting glueball masses are given in Table 3. 
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Table 3 


Glueball masses (in GeV) for L >2, from the Hamiltonian (2) 


L,5 

States 

Masses 

2,2 

0++,4++,l++,2++,3++ 

3.11-3.17 

2,0 

2++ 

3.13 

3,1 

4“+,3"+,2“+ 

3.53 

4,2 

6++ - 2++ 

3.88 

4,0 

4++ 

3.88 

5.1 

6“+,5"+,4“+ 

4.21 


These masses can be compared to (scarce) lattice values. E.g. = 

3.28 ±0.2 GeV in [221, 3.81±0.31 in |2SI and 3.26±0.21 in |21; Mz„t(4++) = 
3.65 GeV in [2] and around 3.5 GeV in [T]; M;at(6++) = 4.2 ± 0.2) GeV in 

Cl 

One can notice an approximate degeneracy of O’*"’'' and states on 
the lattice in agreement with theory. For the pomeron trajectory the state 
6'''+(L = 4, S' = 2) is important and will be used below. 

It follows from Table 3, that there are several states with the same 
which do not differ substantially in mass (for example 4'^^ in (2,2), (4,2) and 
(4,0) conhgurations). This makes difficult to compare predicted masses for 
these states with lattice data, which are usually approximated by a single 
state. Note that masses of states with on the pomeron trajectory 

calculated in ref pQ are higher than the lowest values, obtained in our paper. 
As a result the slope of the pomeron trajectory obtained in paper pQ is 
smaller than the value corresponding to Gasimir scaling and an intercept 
of the purely glueball pomeron is higher than in the present approach. As 
we discuss in Section 4 below this intercept has little to do with a physical 
pomeron intercept due to large mixing with qq trajectories. 


3 Three-gluon glueballs 

The 3g glueballs (oddballs) can be of two basic configurations: the A-type 
and the V-type (only the first one was considered in [S]). The corresponding 
wave operators for the A-type are given in Table 10 of [S] and basically 
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correspond to 3 gluons sitting in vertices of a triangle and connected by the 
fundamental strings. Another, not considered in jH] the Y, form is composed 
as 

'^{xi,X2,X3) ~ dabcE^{x)El{x)Ei{x) (3) 

and is the adjoint equivalent of the baryon operator with the replacement 
^ 0/37 —> dabc- Another possible F-type form is made with the operator fabc 
instead of dabc and requires antisymmetric spin-coordinate function. Using 
the charge-conjugation C transformation ^ ~E^^, one can easily under¬ 
stand that the form made of fabc has C = -|-1. Both /, d forms have been 
used in |3j. 

The spin-independent part of Hamiltonian in both cases can be written 
as 

= T^g + “^ + ^,A(ri5 1 ’ 2 , Ts) (4) 

where T^g is the kinetic operator, 

^ P^ + P| Id Id 

Tsg = —h ;—P = Pg = -tt 

2/i i d^ ' i drj 

and the Jacobi coordinates dehned as 


I 


[3 fri + r2 \ ri - r2 


( 6 ) 


Finally the interaction is 


Va = (TfYl lu - Tjl, Vy = (Jadj X! lu - RyI 

i<j 2=1 


(7) 


and Ry is the position of the string junction. We shall be using a/ = 

C2(adj) 9 

We are solving equation = M(/i)\I/ as in jH] using the hyperspher- 

ical approach jSHISS, which yields very good accuracy already in the lowest 
approximation |HH|. Dehning the hyperradius p, and grand 

orbital momentum K, K = L, L + 2, L -|- 4,.., one has the equation 


1 


( 8 ) 







with 


(9) 


Ua,y{p) — CA,YpO', 

= 2.88; Cy = 3.31. (10) 

Stt 

Here one notice that the case of Y form can be obtained from baryonic 
calculations of |S21 by a simple replacement <Jf ^ <Jadj (see j321 for details 
of derivation) and the total mass M{3g) is obtained by minimizing the mass 
M(/i) over the values of /i, 

M(3(7) = minM(/i) = min + £(/r)l . (11) 

^ L 2 J 

It was shown isa, that e{fi) can be found (with one percent accuracy) from 
the minimum of Ua,y{p) ^-t some point p = po- 
In this way one obtains for e{p). 


£a,y{p) 


3 (ca,y)^/V2/3 ^ 

2 pV3 



pl/3 


( 12 ) 


where k = K‘^ + AK+^. 

Minimizing over p in m, one hnds the constituent gluon mass po, 


/2 

/^o = A/a, M{3g) = 6po. (13) 

One can now predict the 3g glueball masses, still without spin splittings, 
hxing the value of = L = 0,1, 2,... 

For the lowest 3g state with iFmin = L = 0, the spin splitting is due to 
the hyperhne interaction and it was calculated in jH], yielding 

AM^s = —- (14) 

D 

This gives AMss{3 ) = 0.28 GeV, AMss{2 ) = 0, AMss{l ) = 
-0.189 GeV. 

For L > 0 splitting is due to tensor and spin-orbit forces, and we assume, 
that the situation is similar to that of baryons, where these forces are known 
to be weak. The situation might however be different for 3g states. 
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We neglect spin splittings for L > 0 and give in Table 4 the calculated 
spin-averaged masses. One should note that the Y-tjpe glueballs are fully 
equivalent to baryons and their masses are obtained by simply multiplying 

baryon masses from [22] by the factor c2Uund) ~ I sophisticated 

conhgurations were also considered in 133 however the resulting masses are 
very close to the corresponding My values and are omitted from the Table 

4. 


Table 4 

Oddball masses (in GeV) for various L = 0,1, ...4 and J < 7 in compar¬ 
ison with lattice results. For L > 0 only spin-averaged values are shown, 
a = 0.18 GeV2. 


L 

jFT 

M^{theor) 
this work 

Myitheor) 
this work 

M{lat) 

[24] 

[1] 

0 

1— 


3.02 

3.32 

3.40±0.21 

3.10 


2— 


3.21 

3.53 

3.56±0.21 

3.55 


3— 


3.49 

3.83 

3.73±0.21 

4.15 

1 

0+-,. 

..3+- 

3.72 

4.09 

3.46±0.2 

~3.2 

2 

5--,. 

..1— 

4.18 

4.59 



4 

7--,. 

..1— 

4.96 

5.25 




One can see from Table 4, that agreement between theoretical and lattice 
values is reasonable and of the same quality, as the agreement between results 
of different lattice groups. Our results are closer to the lattice data of [221121] 
and he below those of PP; the spin splittings in the (1 ,2 ,3 ) triplet 

are 0.47 GeV in our calculation jS) and 0.33± 0.21 in data of |22l l2Tj . which 
again may indicate that effective for spin-spin interaction is = 0.2 
rather than our hxed value = 0.3. 

Two other models have been used in |2| for 1 , 2 , 3 , 5 , 7 states 

with results similar to ours for the hrst triplet of states, but much heaver for 
5 ,7 states. 
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4 Pomeron and odderon trajectories 

Since our results are the same as in [S] for L = 0, 2, and the resulting pomeron 
trajectory is assumed to be the same as in [H] (mixed with /, f trajectories), 
it is worthwhile to compare the slopes ap(0) obtained from the standard 
procedure, a'p{stand) = = 0.393 GeV“^, with the slopes obtained from 

masses M{J = 2++;4++;6++) = (2.3 - 2.21; 3.15; 3.88) GeV, where we give 
the 2^^ mass interval for ag = 0.3 — 0.2. 

One obtains for the neighboring masses 


ap(4++ - 2++) = (0.431 0.397) GeV^ ap(6++ - 4++) = 0.3971 GeVl 


(15) 


Thus we see that the slopes are close to the standard one, and the trajec¬ 
tory is close to the straight line. The value of intercept, on the other hand 
depends crucially on the intersection with /, /' trajectories and as we argued 
in |H], this value is not actually controlled by the masses on the pomeron 
trajectory. This is in contrast to the case of odderon, where no intersection 
with meson trajectories is possible for t > 0 and hence intercept can be es¬ 
timated from the computed above masses, and now we have three states on 
the trajectory and can determine both slope and intercept. 

We dehne two odderon trajectories corresponding to the A-type and Y- 
type conhgurations. 


ttA(t) = a'^ - t + «a(0), ayit) = a'y ■ t + Q!y(0). 


( 16 ) 


It appears that one can hnd odderon trajectory analytically, using explicit 
mass formulas Eqs. (1I2D-(II1. Here spin splittings are neglected, as it usually 
done for meson Regge trajectories [211 • E-G. in the latter case spin dependent 
correlations split the spin-averaged Regge trajectory into the two nonlinear 
trajectories with J = L ± 1. 

In this way one can approximate k in (fTTIjl for K = L, as k = (L-|-2)^ — | 

(L -|- 2)^ and hence 



Taking into account, that Jodd = L+3, one obtains the odderon trajectory 


M2(3^) = 20.0a(J-aA(0)), aA(0) = 1 - ^ = 0.134. (18) 
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At this point one can compare oddball masses from (jl 8 |l for J = 3 ,5 

and 7 with the valnes given in the Table 4, and hnd that these masses 
differ less than 1%, if for M(3 ) one takes the spin-averaged value, almost 

coinciding with M(2 ) = 3.21. 

It is interesting to compare the resulting slope in 03 “A = si = '>■277 
GeV ^ with the standard relativistic ’’potential” slope (when string rotation 
correction is not taken into account which is our case),a'^{pot) = ^ = 0.308 
GeV“^. The string correction was computed in [171121 and it decreases the 
glueball masses with L 7 ^ 0 by 4 5% and the asymptotic (L S> 1) slope is 

a'^{string) = = 0.39 GeV“^. If one uses this slope (and keeps the mass 

with L = 0 unchanged), then one has the string-corrected odderon trajectory 

A7y.„,(39) = 2’'<^.(.7 - <‘"»(0)), =-1.05, (19) 

It is clear from that masses of F-type glueballs are higher, 

namely =1.1, and the intercept of Y trajectory is lower. 

The same procedure as for Eq. (HI) leads for the intercept of F-type 
trajectory the value 

^^trmg(o) ^ (^ 20 ) 

Now if one insists on using spin-splitted (nonlinear in general) trajectory 
passing through 3 and 5 states, one obtains extrapolating a^it) to 
f = 0 the intercept cta = — 1 . 6 . 

We see that determination of intercepts of 35 >-trajectories is model depen¬ 
dent, but in both cases, Eq. m and Eq. dH), the odderon intercept is far 
from unity. 

5 Discussion and summary 

As shown in the previous section, the odderon intercept is 0.14 (without 
string corrections) and drops down to -1.05 when the string slope is 
restored. In any case intercepts of these 35 >-trajectories are very far from the 
perturbative intercept a^g ~ 1. Thus the nonperturbative effects strongly 
reduce the intercept of the 35 >-trajectories and there is no ’’odderon” (the 
singularity with negative C = a and q;o( 0) ~ 1) in our approach. 

This difference for the intercepts of 35 >-singularities with results of pertur¬ 
bation theory is of principle importance and takes place for all multigluonic 
states. In the perturbation theory the multigluon singularities in j-plane are 
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above unity m and in principle mix with 25 f-state. Account of nonpertur- 
bative interaction between gluons lead to large masses of multigluon states 
and consequently to low intercepts of such trajectories. Qualitatively it can 
be understood as follows: NP interactions conhne and create effective mass 
for each gluon (calculable in FCM [THl[TT)] i proportional to ^/o^. Therefore Ag 
states are separated from 2g states by an interval ~ 2 GeV and mixing can 
be neglected in the hrst approximation. This is contrast to BFKL method, 
where such interval is absent and in principle all multigluon states should be 
considered simultaneously. 

Comparison of our calculated masses with lattice data in Tables 2,3 shows 
that the ordering of states in mass values is the same. The agreement in 
masses for states with L = 0,1 and = 0 is surprisingly good although 
no htting parameters are used. The (small) discrepancy for the states with 

> 0 was discussed above, and it is planned to improve theoretical accuracy 
for L > 0,nr > 0. For L >2 our spin-averaged values are in a good agreement 

with lattice data of [2211211 and lO-i-20% below data of [T1I21I22]- One should 

take into account at this point, that and states occur from L = 2 
and 4, and L = 4 and 6 respectively, and resulting mixing can shift the 
masses from the one-channel values in lattice computations (as well as in our 
analytic approach). 

For ?)g spin-averaged states in Table 3 the agreement with lattice data 
from [221C] is within 10%. 

Our results for the odderon masses of = 5 , 7 are some 0.5 GeV 

below the results of models H^ff and Hm from [2], however qualitatively the 
odderon intercepts are also low (-0.88 and 0.25 respectively). 

In summary, we have calculated high spin states of both 2g and 3g systems 
and found the corresponding masses, in this way extending our results from 
an earlier paper jS] . We have conhrmed our previous results on the pomeron 
trajectoriy. 

In the case of odderon we have found two trajectories, of A and Y type 
with masses differing by ~ 10% and similar intercepts below 0.14. This result 
might explain why odderon is not seen in photonucleon reactions iniEi, and 
is in sharp contrast to the BFKL-type odderon intercept, which is around 1. 
This work is supported by the Federal Program of the Russian Ministry of 
industry. Science and Technology No.40.052.1.1.1112, and by the grant for 
scientihc schools NS-1774. 2003. 2. One of authors (A.B.K.) acknowledges 
partial support of the grants GRDF RUP2-2621-M0-04 and 04-02-17263. 


13 


References 

[1] H.B.Meyer, M.J.Teper, Phys. Lett. B605 (2005) 344. 

[2] Da Quing Liu, Ji Min Wu, hep-lat/0105019, 

[3] F.J.Llanes-Estrada P.Bicudo, S.R.Cotanch, hep-ph/0507205, 

[4] A. Donnachie, H.G.Dosch, O.Nachtmann, hep-ph/0508196, 

[5] C.Ewerz, hep-ph/0511196; hep-ph/0306137, 

[6] C.Adloff et al. HI Collaboration, Phys. Lett. B544 (2002) 35. 

[7] T.Berndt (for the HI Collaboration), Acta Phys. Polonica B33 (2002) 
3499 . 

[8] A.B.Kaidalov, Yu.A.Simonov, Phys. Lett. B 477 (2000) 163; 
hep-ph/9912434; Phys. At. Nucl. 63 (2000) 1428, hep-ph/9911291 

[9] H.Fritzsch and P.Minkowski, Nuovo Cim. A30 (1975) 393; M.Fritzsch 
and M.Cell-Mann, Proc. 16th Int. Conf. on High-Energy Physics, FNAL 
1972 V.2. P. 135. 

[10] R.L.Jaffe and K.Johnson, Phys. Lett. B60 (1976) 201. 

[11] P.Hasenfrantz and J.Kuti, Phys. Rep. 40 (1978) 75. 

[12] D.Robson Nucl. Phys. B130 (1977) 328; 

J.M.Cornwall and A.Soni, Phys. Lett. B120 (1983) 431. 

[13] F.E.Low, Phys. Rev. D12 (1975) 163; 

S.Nussinov, Phys. Rev. Lett. 34 (1975) 1286. 

[14] V.S.Fadin, E.A.Kuraev, L.N.Lipatov, Sov. Phys. JETP 44 (1976) 443; 
45 (1977) 199; LI.Balitsky, L.N.Lipatov, Sov. J. Nucl. Phys. 28 (1978) 
822; L.N.Lipatov, Nucl. Phys. B365 (1991) 614; Sov. Phys. JETP 
63 (1986) 904; J.Bartels, Nucl. Phys. B175 (1980) 365; J.Kwiecinski, 
M.Praszalowicz, Phys. Lett. B94 (1980) 413. 

[15] H.C.Dosch, Phys. Lett. B 190 (1987) 177; 

H.C.Dosch and Yu.A.Simonov, Phys. Lett. B 205 (1988) 339; 

Yu.A.Simonov, Nucl. Phys. B 307 (1988) 512. 


14 





[16] A.Di Giacomo, H.G.Dosch, V.I.Shevchenko, Yu.A.Simonov, Phys. Kept. 
372 (2002) 319; hep-ph/0007223 

[17] A.Yu.Dubin, A.B.Kaidalov, and Yu.A.Simonov, Phys. Lett. B323 
(1994) 41; Yad. Fiz. 56 (1993) 213. 

[18] Dan La Gourse and M.G.Olsson, Phys. Rev. D39 (1989) 2751; 
M.G.Olsson, Nuovo Gim., 107 A (1994) 2541. 

[19] Yu. A.Simonov, Phys. Lett. B 249 (1990) 514. 

[20] L.Lukaszuk, B.Nicolescu, Lett. Nuovo. Gim 8 (1973) 405. 

[21] B.Nicolescu, Nucl. Phys. Proc. Suppl. 25B (1992) 142; hep-ph/9911334 

[22] G.Morningstar, M.Peardon, Nucl. Phys. B63 (1998) 22; (Proc. Suppl.); 
Phys. Rev. D60 (1999) 034509. 

[23] M. Teper, hep-th/98121 

[24] Y.Ghen et al., hep-lat/0510074 

[25] F.Brau, G.Semay, hep-ph/0412173 

[26] V.Mathieu, G.Semay, F.Brau, hep-ph/0511210, 

[27] A.P.Szczepaniak, E.S.Swanson, Phys. Lett. B577 (2003) 61. 

[28] Yu.A.Simonov, in; ”QGD; Perturbative or nonperturbative?” 
L.S.Ferreira, P.Noguera, and J.I.Silva-Marcos eds. Proceedings of the 
XVII Quantum School, Lisbon, 29 Sept. - 4 Oct., 1999, p. 60 (World 
Sci., 2000). 

[29] A.M.Badalian, B.L.G. Bakker, Phys. Rev. D66 (2002) 034025, 
hep-ph/0202246 

[30] V.L.Morgunov, A.V.Nefediev, and Yu.A.Simonov, Phys. Lett. B459 
(1999) 653; hep-ph/9906318 

[31] Yu.A.Simonov, Yad.Fiz., 3 (1966) 630; (Sov.J. Nucl. Phys., 3 (1966) 461; 
A.M.Badalian and Yu.A.Simonov, Yad.Fiz. 3 (1966) 1032; 5 (1967) 88; 
(Sov. J. Nucl. Phys. 3 (1966) 755; 5 (1967) 60; F.Galogero and Yu.A. 
Simonov, Phys. Rev., 169 (1968) 789. 


15 




[32] Yu.A.Simonov, Phys. Lett. B228 (1989) 413; M.Fabre de la Ripelle and 
Yu. A.Simonov, Ann. Phys. (N.Y.) 212 (1991) 235. 

[33] B.O.Kerbikov, M.I.Polikarpov, L.V.Shevchenko, Nucl. Phys. B331 
(1990) 19; 

I.M.Narodetskii, M.A.Trusov, Phys. Atom. Nucl. 67 (2004) 762. 

[34] H.J. de Vega, L.N.Lipatov, Phys. Rev. D66 (2002) 074013; 
G.P.Korchemsky, J.Kotanski, A.N.Manashov, Phys. Lett. B583 
(2004)121; Phys. Rev. Lett. 88 (2002) 122002. 


16 



